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O ■ ABSTRACT 
<N . 

Context. The Galactic microquasar GRS 1915-1-105 exhibits at least sixteen types of variability classes. Transitions from one class to 
^ '■ another could take place in a matter of hours. In some of the classes, the spectral state transitions (burst-off to burst-on and vice versa) 

' were found to take place in a matter of few to few tens of seconds. 

Aims. In the literature, there is no attempt to understand in which order these classes were exhibited. Since the observation was not 
CO , continuous, the appearances of these classes seem to be in random order. Our goal is to find a natural sequence of these classes 

and compare with the existing observations. We also wish to present a physical interpretation of the sequence so obtained using two 
component advective flow model of black hole accretion. 

Methods. In the present paper, we compute the ratios of the power-law photons and the black body photons in the spectrum of each 
class and call these ratios as the 'Comptonizing efficiency' (CE). We sequence the classes from the low to the high value of CE. 
. The number of photons were obtained by fitting the spectra of two independent sets of data of each class with disk blackbody and 

i-^ ' power-law components, after making suitable correction for the absorption in the intervening medium. 

Results. We clearly find that each variability class could be characterized by a unique average Comptonizing efficiency. The sequence 
of the classes based on this parameter seem to be corroborated by a handful of the observed transitions caught by Rossi X-ray timing 
?^ , explorer and the Indian payload Indian X-ray Astronomy Experiment and we believe that future observation of the object would show 

^ ■ that the transitions can only take place between consecutive classes in this sequence. Since the power-law photons are produced by 

^ ' inverse Comptonization of the intercepted soft-photons from the Keplerian disk, a change in CE actually corresponds to a change in 

geometry of the Compton cloud. Thus we claim that the size of the Compton cloud gradually rises from very soft class to the very 
hard class. 

(N 
> 

m 

1. Introduction photons, presumably emitted from the post-shock region of the 

'""^ low angular momentum (sub-Keplerian) flow (Chakrabarti & 

^ The enigmatic steflar mass black hole binai-y GRS 1915+105 Manickam, 2000, hereafter CMOO; Rao et al., 2000). 

S ' ^u^'l^rl^rt ^?^^'^ '^''T?nnl'" ^'^'^^ Small scale variabilities of GRS 1915+105 are identified 

^-H the WATCH detectors (Castro-Tirado et al. 1992) as a tran- -^u i i ■ e ■ j- i /m j- » i oaaa 

I . .,..). , with local variation or the inner disk (Nandi et al., 2000, 

sient source with a significant variability in X-ray photon counts ^, , u e n/r ■ i oaaa r • p d n • oaatn 

- - ,^ ^. , , T ■ n^rrnr- , c- , Ac- Chakrabarti & Manickam, 2000; Mighari & Belloni, 2003). 

(Castro-Tirado et al. 1994). In the RXTE era, GRS 1915+105 c i u u . a .1 . .u- u- . u u . a 

^ . , , , . . , , , , , Several observers have reported that this obiect exhibited many 
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was monitored thousands of times in the X-ray band and the . f . , /v j . i innn n v a q 

^ , . ^ . . types or variability classes (Yadav et al., 1999; Rao, Yadav & 

scientific results reveal a unique nature of this compact object, p^^, 2OOO; Belloni et al. 2000, Chakrabarti & Nandi, 2000; 



H ' The radio observation with VLA suggests apparent superlumi- 1' . 1' oaao n j- ' *u f ' u » ' 

rrt ' , ^ . ... „ ,. , . . ^ , . Naik et al. 2002a). Depending on the variation or photon 

nal nature or its radio lets. Radio observation constrains that its . . » u u a /u a n j 

. ■' , ,, ,. counts in ditterent arbitrary energy bands (hardness ratio) and 

maximum distance is no more than 13.5 kpc and that the let 1 1 j- e r^no im c , iac »u v • u-i 

^ . , , ,. , „ color-color diagram of GRS 1915+105, the X-ray variabil- 

axis makes an angle or 70 with the line or sight (Mirabel & ^ e a ^ u ce^ u ■^ a 

_ , . , . ^ b \ ity ot the source was round to have fifteen arbitrarily named 

Rodriguez, 1994). / o s a. 1 o\ ^ i 

(ff, y6, J, 0, 0, X\^X2^Xi^X'i^ 'I' P' ff) classes. In a 

Continuous X-ray observation of GRS 1915+105 reveals 1999 observation of RXTE, the existence of another class w was 

that the X-ray intensity of the source changes peculiarly in a va- reported (Klein-Wolt et al., 2002; Naik et al. 2002a). In the so- 

riety of timescales ranging from seconds to days (Greiner et al., called x (i-e., Xi to Xa) class, the strong variability as is found 

1996, Morgan et al., 1997). Quasi-Periodic Oscillations (QPOs) in other classes is absent. The classes named ;^i,;t'3,y6 and are 

are observed in a wide range of frequencies. QPOs in this source associated with the presence of strong radio jets (Naik & Rao, 

are associated with different types of X-ray variabilities and their 2000, Vadawale et al., 2003). To understand the above features 

timing properties are correlated with spectral features (Muno et from a dynamical point of view, we carried out a correlation 

al., 1999, Sobczak et al., 1999, Rodriguez et al., 2002, Vignarca study in between temporal and spectral features of this source in 

et al., 2003,). The origin of QPO frequencies between 0.5 to 10 different classes. A preliminary report is presented in Pal, Nandi 

Hz is identified to be due to the oscillation of the Comptonized & Chakrabarti (2008, hereafter PNC08). 
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While a large number of papers have been published in the 
literature on GRS 1915+105, to our knowledge, there is no work 
which actually asked the question: are these classes arbitrary, or 
they appear in a given sequence? The problem lies in the fact that 
no satellite continuously observed GRS 1 9 1 5 + 1 05 . Sporadic ob- 
servations caught the object in sporadic classes. In the present 
paper, we try to show that these classes could be parameterized 
by a common parameter, namely, the ratio between the num- 
ber of photons in the power-law component and the black body 
component. We call this as the Comptonizing Efficiency or CE. 
Since the number of photons in the power-law component de- 
pends on the degree of interception by the so-called hot electron 
cloud or Compton cloud (Sunyaev & Titarchuk, 1980, 1985), 
different classes are therefore parameterized by the average size 
of the Compton cloud. Along with the dynamical evolution of 
CE, we compute the spectrum and the power density spectrum 
(PDS) for each of the classes. To accomplish the computation 
of CE, we separate out the photons ysB of the black body com- 
ponent and the photons ypi from the power-law component and 
take the running ratio CE as a function of time to study how the 
Compton cloud itself varies in a short time scale. Our findings 
reveal that the Compton cloud is highly dynamic. We present 
possible scenarios of what might be occurring to it in different 
variability classes. 

The paper is organized as follows: in the next Section, we 
present a general discussion of the observation, our criteria of 
selection of data for analysis and analysis technique. In §3, we 
discuss the procedure of calculation that is adopted to calculate 
the photon numbers. In §4, the results are presented. In §5, we 
present a unifying view where we show that the Comptonizing 
efficiency may be a key factor to distinguish among various 
classes. Finally, in §6, we make concluding remarks. 



2. Observation & Data Analysis 

The RXTE science data is taken from the NASA HEASARC 
data archive for analysis. We have chosen the data procured in 
1996-97 by RXTE as in this period GRS 1915+105 has shown 
almost all types of variabilities in X-rays. Subsequently, in 1999, 
another class, namely o) was seen. In the present paper, since 
we are interested in sequencing these sixteen classes, we will 
rename them as follows: 1=0; 11=6; 111=7; IV=w; V=ju; VI=v; 
Va=A; VIII=/c; IX=p; X=^; XI=a; XII=6i; XIII=;t'2; XIV=;t'4; 
XV=;^'i; XYl=X3). This would be our final sequence also, and 
hence it is easier to remember. Moreover, theoretical discussions 
often uses these Greek characters it would be confusing to use 
these symbols. 

During the data analysis we excluded the data collected for 
elevation angles less than 10°, offset greater than 0.02° and dur- 
ing the South Atlantic Anomaly (SAA) passage. In Table. 1, the 
details of the data selection and ObsIDs are given which we an- 
alyzed in this paper. 

Here we discuss how we analyzed the temporal and spectral 
properties of the data. 

2.1. Timing Analysis 

In timing analysis of the RXTE/PCA data we use "binned mode" 
data which was available for 0-35 channels only, with a time res- 
olution of 2"^ sec and 'event mode' data with a time resolution 
of 2~^ sec for the rest of the channels. We restrict ourselves in the 
energy range of ~ 2 -40 keV for the timing analysis of PCA data. 
After extracting the light curves (using standard tasks) from the 
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Table 1. The ObsIDs and the dates of RXTE data analyzed in 
this paper. * represents the result shown in Fig. 3 of this paper. 



two sets of different modes, we add them by the FTOOLS task 
"Icmath" to have the whole energy range light curve (2 to 40 
keV). The Power Density Spectrum (PDS) is generated by the 
standard FTOOLS task "powspec" with normalization. Data is 
re-binned to 0.01s time resolution to obtain a Nyquist frequency 
of 50 Hz as the power beyond this is found to be insignificant. 
PDSs are normalized to give the squared rms fractional variabil- 
ity per Hertz. 

To have the timing evolution of QPO features in different 
classes, we generate PCA light curves of 2 to 40keV photons 
with 0.01s time bin. We use the science data accordingly from 
the 'binned mode' and the 'event mode' for 10s time interval 
in each step. This light curve is then used to make PDS using 
the powspec task. The dynamic PDS is made by considering the 
shift of time interval by 1 s. The selection of time interval is done 
by FTOOLS task 'timetrans'. The PDSs are plotted accordingly 
to see the variation of QPO features in a particular class of a few 
hundred seconds of observation. 

2.2. Spectral Analysis 

Spectral analysis for the PCA data is done by using "standard2" 
mode data which have 16 sec time resolution and we constrained 
our energy selection up to 40 keV to match with the timing 
analysis. The source spectrum is generated using FTOOLS task 
"SAEXTRCT" with 16 sec time bin from "staiidard2" data. The 
background fits file is generated from the "standard2" fits file by 
the FTOOLS task "runpcabackest" with the standard FILTER 
file provided with the package. The background source spec- 
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tram is generated using FTOOLS task "SAEXTRCT" with 16 
sec time bin from background fits file. The standard FTOOLS 
task "pcarsp" is used to generate the response file with appro- 
priate detector information. The spectral analysis and modeling 
was performed using XSPEC (v. 12) astrophysical fitting pack- 
age. For the model fitting of PCA spectra, we have used a sys- 
tematic error of 1%. The spectra are fitted with diskbb and power 
law model along with 6.0 x lO^^cm"^ hydrogen column absorp- 
tion (Muno et al., 1999) and we used the Gaussian for iron line 
as required for best fitting. During fitting of the spectra we have 
adopted the method used by Sobczak et al., (1999), to obtain 
the values of spectral parameters. We have calculated error-bars 
of 90% confidence level in each case. In case of soft states, the 
diskbb is found be fitted within 2.5 - 10.0 keV. But in case of 
hard states, the diskbb is fitted within 2.5 - 5.0 keV. This up- 
per limit changes dynamically for the fitting of spectra. Table. 2 
gives the energy up to which the black body fit is made for each 
class. 

To have the spectral evolution with time for each class, we 
have generated PCA spectram (2.5 to 40keV) with a minimum of 
16s time interval along with background spectram and response 
matrix. This procedure is repeated with every 16s shift in time 
interval since the minimum time resolution in 'standard2' data 
is 16s. We use 'timetrans' task to select each and every step of 
time interval. Spectral evolution (dynamic representation) for a 
longer duration of observation (~ 500 to 2000 sec) of each class 
is plotted to show the local spectral variation in each class. 

3. Calculation Procedure 

In this section, we discuss the process of calculation of photon 
numbers using the parameters obtained from spectral fit. 

Each spectram is first fitted within XSPEC envirormient with 
the models and constraints mentioned above. The fitting param- 
eters are used to calculate the black body photon counts from 
the Keplerian component and the power-law photon counts from 
the Compton cloud. The number of black body photons are ob- 
tained from the fitted parameters of the multi-color disk black 
body model (Makishima et al., 1986). This is given by, 



f{E) 



—rl cos i {TITi„r''i'B{E, T) dT/Ti„, 



where, B(E, T) - ^^^^g^^ ^^ and r,„ can be calculated from, 

K = in„/iD/lOkpc)fcosi, 



(1) 



(2) 



where, K is the normaUzation of the blackbody spectram ob- 
tained after fitting, r,„ is the inner radius of the accretion disk in 
km. Tin is the temperature at r,„ in keV. D is the distance of the 
compact object in kpc and ; is the inclination angle of the accre- 
tion disk. Here, both the energy and the temperature are in keV. 
The black body flux, /(£) in counts/s/keV is integrated between 
0.1 keV to the maximum energy (obtained from diskbb model 
fits with the spectra as done in Sobczak et al., 1999). This gives 
us with ybb, file number of the black body photons in counts/s. 

The Comptonized photons jpl tiiat are produced due to 
inverse-Comptonization of the soft black body photons by 'hot' 
electrons in the Compton cloud are calculated by fitting with the 
power-law given below. 



PiE) = NiE)-", 



(3) 



where, a is the power-law index and A' is the total 
counts/s/keV at IkeV. It is reported in Titarchuk (1994), that the 



Comptonization spectram will have a peak at around 3 X r,„. 
The power law equation is integrated from 3 x r,„ to 40keV to 
calculate total number of Comptonized photons in counts/s. The 
Gaussian function Ga(E) incorporated with power law equation 
for the presence of line emissions in the spectrum, whenever nec- 
essary is given below. 



Ga{E) = L — exp(-0.5 {e - E,)lcrf)\ , 



(4) 



where, Ei is the energy of emitted line in keV, <j is the line width 
in keV and L is the total number of counts/s in the line. The vari- 
ation of the Comptonizing Efficiency (CE) with time is plotted 
to have an idea of how the Compton cloud may be changing its 
geometry. 

We fitted the data with another Comptonization model, 
namely, the compst (Sunyaev & Titarchuk, 1980). The spectram 
is given by. 



Cst{x) = 



a+3 



aja + 3) / X u ^ 

— ' when < X < xo, 

2a + 3 \ xq / 

a(a + 3)xo" 



r(2Qr -I- 4) 



x^ exp(-x)/(x), when x> xq. 



(5) 



J "CO 1 / , \ . - _ 

Q f""' exp-f(l -H ^] , X = E/kTe, xo 

kTin/kTe, and a is the spectral index. This a is calculated from 

- yjj + y ~ I' where y = ^^r^j^- is the electron tem- 
perature and T is the optical depth of the electron cloud. In the 
computation, we exclude the hydrogen column absorption fea- 
ture while calculating CE. This is because we are interested in 
the photons which were emitted from the disk before suffering 
any absorption. Thus the variation of CE along with other fea- 
tures (e.g., spectral and QPO frequency variations) will reflect 
the actual radiative properties of the flow near a black hole. 

In Fig. 1 , we show an example of the fitted PCA and HEXTE 
spectram along with the fitted components. We also show the 
residuals to characterize the goodness of the fit. In the left panel, 
the fitting is done with diskbb and power-law components. In the 
right panel, the same spectram is fitted with diskbb and compST 
models. The photon numbers and CE are calculated with the 
parameters obtained from the fitting. The black body spectrum 
is simulated with r,„ = 1.15 1q q52 keV. The calculated num- 
ber of black body photons from 0.1 - 6.5keV is 156.03 
kcounts/s. In the same way, the power-law spectram is simulated 
with the power-law index = 2.057 [47- The calculated number 
of the Comptonized photons between 3.45 - 40 keV is 0.39 +^ ^5 
kcounts/s. The ratio between the power-law photon and the black 
body photon is 0.25 "7%. This means that only 0.41% of the 
soft photons are Comptonized by the Compton cloud. For the 
sake of comparison, the simulation with the compST model pa- 
rameters are found to be KTe =6.41 +"5 and t = 9.42 '^l-^. 
The number of Comptonized photons wiihin the same range as 
previous case is 0.41 ^^ '^^ kcounts/s. Thus the ratio between the 
C0OT/75r photon and the diskbb photon is 0.26 ;!;q q7%. In Fig. la, 
we also include HEXTE data, and fitted with a power-law. But 
the spectral fit does not change in slope and the number of pho- 
tons contributed by HEXTE regime is so low that the basic re- 
sult of CE is not changed. In Fig. lb, we give an example where 
comptST model cannot be fitted with HEXTE data. This is true 
for many variability classes. We therefore ignore the HEXTE 
data in the rest of the paper. 

In Table. 2, a comparison of the fitted parameters for all the 
classes is provided with error at 90% confidence level. The first 
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column gives the class name and the second column gives the 
burst-off (harder) and the burst-on (softer) states, if present, in 
each class. These states are actually similar to hard and soft 
states of a black hole, except that the wind plays a major role 
and thus the states are not all the way hard or soft. The third 
column gives the black body temperature of the fitted Keplerian 
disk (Tin in keV). The fourth column gives the maximum energy 
in keV upto which diskbb model is fitted and soft photon num- 
ber is calculated. The fifth column gives the count rates of the 
soft photons. The next column gives the index associated with 
the fitted power-law component. The seventh column gives the 
derived hard photon counts per second. The eighth column gives 
the ratio of the photon numbers given in the fifth and the seventh 
columns. This is the so-called Comptonizing efficiency or CE. 
We arranged the classes in a way that the average CE changes 
monotonically. We find it to be high in harder classes such as 
XIII-XVI and low in the softer classes such as I-II. The final 
column gives the value of reduced for the fits. The behavior 
of average CE will be discussed later. 

As mentioned earlier, we fit the spectrum of each class with 
diskbb and compST model to have the information of r and 
of electron distribution inside the Compton cloud. One sample 
spectrum is taken for each class and fitted with the diskbb and 
compST model with hh — 6.0 x \(fi^crrr^ hydrogen column for 
absorption and 1 % systematic error. The best fitted parameters 
are given in Table. 3 for comparison with the parameters given 
in Table. 2. We find that in the cases where both models were 
used, the CE obtained were very close to each other. 

4. Results 

We discuss below the main results for each class one by one. The 
Observational IDs are given in Table 1 . So far, we have shown 
that the CE, which were computed model independent way, vary 
from one variability class to another, i.e., our results did not de- 
pend on any specific theoretical model. However, to interpret the 
results, and to facilitate the discussions, it is instructive to keep a 
paradigm in mind. In Fig. 2(a-c), we present a cartoon diagram 
which is basically the two component advective flow (TCAF) 
model of Chakrabarti & Titarchuk (1995, hereafter CT95), suit- 
ably modified to include outflows as in Chakrabarti & Nandi 
(2000) (Inclusion of the jets and outflows does not make our 
model a three component model, since the outflows are results of 
the inflow and formed from the post-shock region, which is the 
CENtrifugal pressure supported BOundary Layer, or CENBOL.) 
In the Figure, the BH represents the black hole, the dark shaded 
disk is the Keplerian flow and the light shaded disk is the sub- 



Keplerian flow. The CENBOL and the outflow from it intercept 
the soft photons from the Keplerian disks and may be cooled 
down if the Keplerian rate is sufficient. 

In the cartoon diagram of Fig. 2a, the CENBOL is not cooled 
enough and the Comptonization of the soft photon is done both 
by the CENBOL as well as the outflow. This configuration 
typically produces a 'harder' state. However, the outflows de- 
pends on the shock strength (Chakrabarti, 1999, hereafter C99) 
and could produce weak jets as in XIII-XIV or strong jets as 
in XV and XVI (see, Vadawale et al., 2001). In Fig. 2b, the 
Comptonization is due to the high accretion rate in the Keplerian 
disk and the CENBOL collapses. As a result, no jets or out- 
flows are formed. This configuration typically produces the soft 
states. If both the Keplerian and the sub-Keplerian rates are com- 
parable, then the intermediate states would be produced also. 
In Fig. 2c, we show the situation where the shock strength is 
intermediate and consequently, the outflow rate is the highest 
(see, C99; Chakrabarti & Nandi, 2000 and references therein). 
In this case, there is a possibility that the outflow may be cooled 
down by Comptonization if the intercepted soft photon flux is 
high enough and the outflow is temporarily terminated. The flow 
falls back to the accretion disk, increasing local accretion rate 
in a very short time scale (seconds). We believe that the vari- 
abiUty classes (such as V-XII) some of which show clear softer 
(burst-on) and harder (burst-off) states alternately and showed 
evidences of intermittent jets (Klein- Wolt et al. 2002, Rodriguez 
et al. 2008) belong to this category (Chakrabarti & Manickam, 
2000). 

We now present the dynamical analysis of the light curves of 
all the variability classes. While choosing the data duration for 
a given class, the following considerations have been followed: 
If the count rate has no obvious rise/fall signatures, we use the 
data of 500s. If the count rates have some 'repetitive' behavior, 
we use the data of 500, 1000, 1500s or even 2000s, which ever 
is bigger so that at least one full cycle in included in the data 
chunk. However, in the latter cases, while computing the CE, we 
use the data from a full cycle only. 

In each of Fig. 3(a-p) we show four panels. The top panel is 
the variation of the photon counts with time. The second panel 
is the variation of log(power) of the dynamical Power Density 
Spectra (PDS) which may show presence or absence of quasi- 
periodic variations or QPOs. The third panel is the dynamic en- 
ergy spectrum which shows whether the variability class is dom- 
inated by hard photons or soft photons. Finally, the bottom panel 
shows the variation of Comptonizing Efficiency CE calculated 
using 16 seconds of birmed data. The error bar is provided at 
90% confidence level. 
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Table 2. Parameters for the simulated spectral fits with diskbb plus powerlaw models for aU the variability classes. 
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Table 3. Parameters for the simulated spectral fits with diskbb and compr^r model for all the variability classes. Those classes which 
could not be fitted with this model have been marked with dashes 



4.1. Class No. I 

The Class 1 data shows a very short time scale variability with 
the presence of short scale dips in its hght curve. The results are 
shown in Fig. 3a. Dynamic PDS shows no signature of QPOs 
and the spectrum is mostly soft in nature. The soft photon num- 
ber varies at around 330.24 kcounts/sec while the Comptonized 
photon number varies around 0.16 kcounts/sec. The CE is only 



around 0.05%. Thus in this class very few black body photons 
are intercepted by the CENBOL and correspond to a situation 
similar to that in Fig. 2b. 

4.2. Class No. II 

In Fig. 3b, we show the result of our analysis of the data that 
belong to the Class II. In this case, there is an absence of QPO 
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Fig.2. Cartoon diagrams for three major types of variabilities in GRS 1915+105. The classes XIII-XVI belong to the group (a) 
where the Keplerian rate is low, the CENBOL is large and < CE > is high. Jet/outflow rates may be low and continuous. This is 
the so-called 'Hard' group. The softer classes such as I-IV belong to group (b) where the CENBOL is very small and very little 
Comptonization may be due to the sub-Keplerian flow only. This is the 'Soft' group. The rest of the classes belong to the variations 
of group (c) or the 'Intermediate' group, where, the jet also plays dynamically important role in shaping the spectrum. 



in the dynamic PDS and the spectrum is soft, mostly dominated 
by the blackbody photons. The CE is around 0.06%, indicating 
a situation similar to Class I. 



is not really formed, i.e., the sub-Keplerian flow has very low 
energy and/or angular momentum (C90). 



4.3. Class III 

The class III data appears to be less variable with a distinct and 
repeated 'dip' like features at a time gap of a few seconds. The 
result is shown in Fig. 3c. The class is a steady soft state with no 
QPO visible in the PDS. The CE varies around 0.06% indicating 
the low interception of the soft photons with the CENBOL. 



4.6. Class VI 

The analysis of 500s of class VI data is shown in Fig. 3f. The 
fluctuations of the photon count rate and the spectrum are erratic. 
QPOs are visible and CE is higher only when the spectrum is 
harder. CE varies between 0.04% to ~ 0.3%. 



4.7. Class VII 



4.4. Class IV 

The results of the analysis of the data of 1000s in class IV is 
shown in Fig. 3d. The photon count seems to be steady at around 
10000 counts/sec for most of the time but sometimes for a du- 
ration of around 100 sec the photon count is decreased to about 
2000 counts/sec. In the whole class, no QPO is observed and 
whenever the photon count is low, the spectrum appears to be- 
come harder. The CE remains low at around 0.07. 

The four classes I-FV belong to the softer class. The CE is 
very low, even when the Keplerian rate is high (spectrum is dom- 
inated by soft photons). This means that the CENBOL is very 
small in size and hence the QPOs are also absent. In our picture, 
these cases would correspond to that in Fig. 2b, only accretion 
rates vary. From the count rates, it appears that the disk rates 
are intermediate in classes 1 and II while it is high in class III. 
The sub-Keplerian rate is low in class II, but is intermediate in 
classes I and III. 



4.5. Class V 

The result of our analysis for 500s of data of class V is shown in 
Fig. 3e. In this class, the spectrum, while remaining soft shows 
a considerable fluctuation. CE remains low at around 0.05% to 
0.15%. It goes up when the spectrum is harder. This shows that 
the Keplerian rate may be changing rapidly and the CENBOL 



The result of our analysis of 1000 seconds of data of class Vll 
is shown in Fig. 3g. For the first 200 - 600 sec the photon num- 
ber varies between 10000 and 30000 counts/sec. The black body 
photon rate of the simulated spectrum is around 300 kcounts/sec 
and Comptonized photon rate is around 0.19 kcounts/sec. This 
class is a mixture of the burst-off and burst-on states. When the 
photon count is higher, the spectrum is softer and the object is in 
the burst-on state and when the photon count rate is lower, the 
spectrum is harder and object is in burst-off state. The CE factor 
varies between 0.05% to ~ 0.1%. Here too, the physics of vary- 
ing the Comptonizing region (and thus the CE factor) is similar 
to what is seen in classes VIII-IX below. At the end of the 400s 
span when the burst-off state starts, a strong spike in both the 
CE factor and QPO frequency are observed, indicating a sudden 
flaring in the CENBOL configuration. 

4.8. Class VIII 

An analysis of a 500s chunk of the data of class VIII of obser- 
vation is shown in Fig. 3h. In the class VIII, the photon counts 
become high ~ 30000/sec and low (~ 10000/sec) aperiodically 
at an interval of about 50 - 75s. In the low count regions, the 
spectrum is harder and the object is in the burst-off state. Distinct 
QPOs are present and at the same time, the Comptonizing effi- 
ciency is intermediate, being neither as high as in the class no. 
XIII-XVI, nor as low as in class 1-111. In CENBOL picture, low 
frequency OPOs are indicators of the shock oscillations (CMOO). 
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Fig. 3. Results of the analysis of (a) I, (b) II, (c) III, (d) IV,(e) V, (f) VI, (g) VII, (h) VIII, (i) IX, Q) X, (k) XI, (1) XII, (m) XIII, (n) 
XIV, (o) XV, (p) XVI classes are shown. First Panel: Light curve in 2.0-40.0 keV range. Middle Panel: Dynamic PDS of the light 
curve. Clear evidence of low frequency noise and QPOs are seen. Third panel: Dynamic PCA spectrum showing subtle variations 
of the spectral characteristics with time. Bottom Panel: Comptonizing Efficiency (CE) in percent, obtained from 16s binned data. 



As discussed before, in this case the strength of the shock is 
intermediate and produces strong outflows (CMOO). In the fit- 
ted spectrum, the Keplerian photon varies between 150 to 450 
kcounts/sec and Comptonized photon varies between 0.39 to 
0.21 kcounts/sec. The CE factor becomes high at 0.25%. The 
CENBOL is cooled down. After the matter which fell back on 
the CENBOL is totally drained, the burst-ofl^ state is resurfaced 
with a lower CE (less than 0.1%). 

4.9. IX class 

The result of our analysis of the class IX data is shown in Fig. 3i. 
This class contains a cychc variation of the hard and soft photons 
with roughly 75 s of periodicity. The blackbody photon count in 
the simulated spectrum varies between 350 to 420 kcounts/s and 
Comptonized photon counts vary around 0.95 to 0.34 kcounts/s. 
In the harder state (low count regions), the CENBOL is promi- 
nent and the Keplerian component is farther away. Thus the QPO 
is prominent there also. A gradual softening of the hard spectrum 
indicates that the Keplerian disk is moving towards the black 
hole and the CENBOL is becoming smaller in size while still re- 



taining its identity. Keplerian photon interaction increases with 
CENBOL size and CE rises to a maximum of 0.25%. At the peak 
of the light curve, the CENBOL, which is also the base of the jet 
is cooled down due to the increased optical depth. 



4.10. Class X 

The result of analysis of 1600s of Class X data is shown in Fig. 
3j. In the first phase of 600 sec, the quasi periodic variation of 
photon counts takes place in gradually decreasing period from 
10000 counts/sec to 30000 counts/sec. In this phase, the spec- 
trum is soft and the QPO is seen only when the photon count is 
low. Here, the CENBOL is visible when the Keplerian flow is 
going far away from time to time. The CE varies around 0.14%. 

In next 800 sec the spectrum is harder and the distinct varia- 
tion of QPO frequency (from 12 Hz to 3 Hz) indicates the varia- 
tion of the shock location. However, in this phase, CE is varying 
around 0.27%. This means that outflow is taking an active role 
in intercepting the soft-photons. 
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4.11. Class XI 

Class XI is an intermediate class in which QPO is always ob- 
served, the QPO frequency is generally correlated with the count 
rate (as in XIII-XVI classes discussed below). Since this class 
displays a long time variation, we analyze the data of 2000s. Fig. 
3k shows the result. The CE varies between 0.05% to ~ 0.6%. 
The spectrum shows that the Keplerian rate is not changing 
much, but the sub-Keplerian flow fluctuates, perhaps due to 
failed attempt to produce sporadic jets. The CE varies consid- 
erably and so does the CENBOL. 



4.12. Class XII 

The result of the analysis of a 1000s data of class XII is shown 
in Fig. 31. This class can be divided in two regions depending on 
the photon count rates. In the soft dip region the photon count 
rate is lower than 10000 counts/sec. In this region, the spectrum 
is softer and the CE amount of interaction is around 0.17%. 

In the other (hard dip) region, say between 350s and 820s, 
the photon count is higher and varies from 10000 counts/sec to 
30000 counts/sec. In this region, CE reaches a high value of 
0.28%. The QPO is also present. Here, the CE factor, which 
is linked to the geometry of the Comptonizing region, is anti- 
correlated with the QPO frequency (and correlated with shock 
location, i.e., CENBOL size). In the third panel, where the dy- 
namic spectrum is drawn, we note that between the soft dip 
and the hard dip, the low energy photon intensity is not chang- 
ing much, while the intensity of spectrum at higher energies is 
higher. This indicates that the sub-Keplerian rate is high in the 
hard dip state. 



4. 13. Classes No. XIII & XIV 

Xin and XrV classes have low soft X-ray fluxes with less intense 
radio emissions. Dynamic PDS shows a distinct QPO at around 3 
Hz in class XIII and around 5 Hz in class XIV. The X-ray photon 
is significant up to 40 keV with a flatter power-law slope, which 
signifies that the source belongs to a hard state. The results are 
given in Fig. 3m and Fig. 3n. The Comptonizing efiiciency in 
XIII and XIV classes is around 0.4 - 0.5%. 



4. 14. Classes no. XV & XVI 

The classes XV and XVI correspond to the radio-loud states, 
whereas classes XIII and XIV are in radio-quiet states (Vadawale 
et. al., 2001). In both the cases, the dynamic PDSs show a strong 
QPO feature around 4 Hz and 0.8 Hz respectively through- 
out the particular observation and the spectrum is dominated 
by the hard power law. Results are shown in Fig. 3o and Fig. 
3p. Around 0.77% of the soft photons are intercepted by the 
CENBOL in XIII class, whereas it is 0.88% in XVI class. In 
our picture, the CENBOL (Fig. 2a) oscillates to produce the ob- 
served QPOs (Chakrabarti, 1997; CMOO). Generally, XIII-XVI 
classes are in harder states and in the language of two compo- 
nent model (CT95, C97), the Keplerian rate is low and unable to 
cool the CENBOL and jet combined which are produced by the 
sub-Keplerian flows. Whether the jets would be produced or not 
depends on the shock strength (C99). 



5. Comptonizing efficiency in different classes: an 
unifying view 

In the above sections we have analyzed the data of all the six- 
teen classes. We presented the light curves, the dynamical en- 
ergy and power density spectra and the ratio of the power-law 
photons to the black body photons from the simulated fit of the 
original spectra. In the literature, there has been so far no discus- 
sion on the sequence in which the class transitions should takes 
place. Also, there is no discussion on which physical properties 
of the flow, the characteristics of the variability classes depend. 
In the present paper, we analyzed many aspects of the variabiUty 
classes which may be understood physically from the two com- 
ponent advective disk paradigm. Since CE behaves diff'erently in 
various classes, it may play an important role in distinguishing 
various classes. This CE is directly related to the size and optical 
depth of the Compton cloud which determines the fraction of the 
injected soft photon which are intercepted by the Comptonizing 
region. 

Though CE varies very much in a given spectral class, it is 
instructive to compute the average CE in a given class (aver- 
aged over a period characterizing the class). We denote this as 
< CE >. In Fig. 4, we present the variation of log(< CE >) in 
various classes. We placed error bars also in all the average val- 
ues which are at 90% confidence level. We arranged arbitrarily 
defined sixteen classes in a manner so that < CE >s are mono- 
tonically increasing. This gives rise to a sequence of the variabil- 
ity class shown in the X-axis. In classes X and XII we did not av- 
erage CE over both hard and soft regions, since it is believed that 
the soft regions are produced due to sudden disappearance of the 
hard region (e.g., Nandi et al. 2001). Thus, while placing them 
in the plot, we considered the average over the burst-off (harder) 
state only. To show that the sequence drawn is unique, we plot- 
ted the average values of two sets of data of aU the classes. The 
average of the first set is drawn with a dark square sign and the 
average of the second set is drawn with a dark triangle sign. The 
individual error bars are drawn with solid and dashed lines re- 
spectively. In both the sets, the sequence is identical with softer 
classes to the left and harder classes to the right. Based on the 
nature of the variations of CE, we divided these classes into three 
groups corresponding to three types of accretion shown in Fig. 
2(a-c). Classes I to IV belong to 'softer' state and classes XIII- 
XVI belong to 'harder' state. The rest of the classes belong to 
the 'intermediate' state. 

The above findings of alignment in average CE variation 
shows that it plays a significant role in distinguishing the vari- 
ability classes. Moreover, as we discussed already, the variation 
in CE for a given class is due to the softening and hardening of 
the spectra in a much shorter time-scale and thus is related to 'lo- 
cal' physical processes, such as interaction with winds and out- 
flows, which, in turn, depends on the existence or non-existence 
of CENBOL, i.e., the spectral states. 

In Fig. 5 (a-b), we plotted the results of these two sets of 
analysis separately. We show separate averages over CE in the 
burst-ofl" and burst-on states (as two dark circles) when they are 
present as well as the global average over CE in a given class 
by filled squares (except classes X and XII where the physics is 
diff'erent and it is meaningless to talk about the overall average; 
and hence only average over the harder region is plotted.). 

In Chakrabarti et al. (2004) and Chakrabarti et al. (2005), it 
was mentioned that although many observations were made of 
GRS 1915-1-105, only a few cases, direct transitions were ob- 
served. In particular, using Indian X-ray Astronomy Experiment 
(IXAE) they showed direct evidences of /c — > p (VIII IX), 
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Fig. 4. Variation of average Comptonizing efficiency (< CE >) for different variability classes of GRS 1915+105. The error-bars 
were calculated from the excursion of CE in a given class. Filled squares and triangles represent < CE > for two sets of variabiUty 
classes. Generally softer classes have smaller < CE > and harder classes have higher < CE >. 
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Fig. 5. Variation of averaged Comptonizing efficiency (< CE >) for different variabihty classes of GRS 1915+105. In the left panel, 
we show the results for the data set plotted with filled squares in Fig. 4, and the right panel, we show the results of the data plotted 
with filled triangles in Fig. 4. The averages in hard and soft chunks are placed separately with filled circles. In class X averages of 
three chunks (burst on, pre-spike and post-spike) are plotted with filled circle, filled box and filled circle respectively. The sequence 
is found to remain the same in both sets of data. We placed classes X and XII according to < CE > in the hard (pre-spike) region 
since the softer regions are beUeved to be produced due to totally different physical processes. 



X ^ P (XIII ^ IX), X (XIII^ Xll) and p ^ a (IX^ XI) 
transitions in a matter of hours. In Naik et al. (2002), IXAE data 
was used to argue that p class variability could have changed 
to X class via a class (i.e., IX ^ XI ^ XII). In Nandi et al. 
(2001), it was shown that the 6 class (Class XII) is rare and the 
observed soft-dip is perhaps due to the disappearance of the in- 
ner region by magnetic rubber band effect. Accepting that class 
xn as anomalous, we find from Fig. 4 that the observed transi- 
tions reported in Chakrabarti et al. (2004, 2005) and Naik et al, 
2002 are 'naturally' explained. For instance, — > p (0 and y6 be- 
ing anomalous and intermediate a has been reported by Naik et 
al. 2002), p ^ a (fi being anomalous) and k ^ p, are expected 
from our analysis. Similarly, we can claim that there should not 
be any transitions such as I ^ IX; III^ VII; for example. We be- 
lieve that if we carry out spectro-photometry of GRS 1 9 1 5 + 1 05 
continuously (Chakrabarti et al. 2008), then we may be able to 
catch the transition from one type to another more often and ver- 
ify if the sequences we mentioned here need further refining. 



6. Conclusions 

In this paper, we have analyzed all known types of light curves 
of the enigmatic black hole GRS 19 15 +105 and computed the 

dynamical nature of the energy and the power density spectra. 
We did not characterize these classes by conventional means. 



such as using hardness ratios defined in certain energy range 
since such a characterization does not improve our view about 
the physical picture. Furthermore, characterization using certain 
energy range is possible only in a case by case basis, and is not 
valid for the black holes of all masses. Instead, we asked our- 
selves whether we can distinguish one class from another from 
physical point of view purely in a model independent way. We 
observed that independent of what the nature of the Compton 
cloud is, the weighted mean < CE > of Comptonizing efficiency 
(CE) obtained every 16 seconds of the binned data, increases 
monotonically as the class varies. This pattern we find was veri- 
fied with two sets of data covering all the variability classes. So, 
iCE^ is not arbitrary - it is characteristics of a class. Based on 
the values of jCE^, it is observed that the classes belong to three 
states: Classes I-IV in softer states. Classes V-XII in intermedi- 
ate states, and Classes XIII-XVI belong to the harder state. 

When the weighted average value of CE is monotonically ar- 
ranged, we obtain a sequence which appears to be followed dur- 
ing the actual transitions. Indeed, when comparing with avail- 
able data of PCU of RXTE and the Indian payload IXAE on- 
board IRS-P3, the transitions from one class to another as re- 
ported in the Uterature do follow our sequence. Given that a large 
variation of CE occurs in a given variability class, it is puzzling 
why the sequence of the classes obtained by us should follow 
changing the CE values averaged over the whole class. It is pos- 
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sible that a specific value of < CE > actually forces the system 
to be in a given class, just as a parameter such as wind speed 
decides the mean angle of an oscillating pendulum. The excur- 
sion of CE in that class could be due to totally different physical 
process and not necessarily due to interception of soft photons 
by the CENBOL alone. May be the mass-loss rate of CENBOL 
is playing a role, which in turn depends on the shock strength. 

In the two component model of the Chakrabarti-Titarchuk 
(CT95), the CENBOL and associated outflow play the role of 
the Compton cloud. In this model, the increase in Compton ef- 
ficiency can be affected in several ways: (i) by increasing the 
shock location increases the size of the CENBOL which inter- 
cepts larger number of soft photons and/or (ii) by increasing 
the accretion rate of the sub-Keplerian component, which in- 
creases the optical depth and scatter more soft photons to pro- 
duce power-law photons. In CT95 and C97 it was shown that 
harder states are produced by both the effects mentioned above. 
Indeed, we find that the classes with harder states have more CE. 
Shifting of the shock locations is possible by changes in viscos- 
ity or changes in cooling rate in the post-shock region. The time 
scales of such effects in a sub-Keplerian flow could take hours. 
Time scale of changing the global sub-Keplerian flow rates could 
be comparable to the free-fall time from the outer edge, i.e., of 
the order of a day or so. On the other hand, in Classes VII and 
VIII etc., the CE changes in a matter of minutes. Such a short 
time variability of Comptonizing efficiency is possible if the base 
of the outflow is abruptly cooled and returned back to the disk in- 
creasing the accretion rate of the Keplerian/sub-Keplerian rates 
locally (CMOO). 

In this paper, we find that the < CE > is really important 
in deciding the sequence. However, for a given class, the degree 
of excursion of CE in a given variability type must depend on 
another parameter, such as the outflow rate which in turn de- 
pends on the shock strength. A simple estimate (C99) suggests 
that the shock strength decides the outflow rate and hence the 
time taken by the base of the outflow to reach unit optical depth 
T for Compton scattering (CMOO). For a very strong shock, this 
outflow rate is very weak, as is evidenced by weak (few tens of 
miliJansky) radio flux even for 'radio-loud' classes. For interme- 
diate shock strength the outflow rate is higher, and it is easier to 
have T = 1 in a short time scale (CMOO) and fractional change 
in CE also becomes high. In classes XII-XVI, we not only see 
CE to be very high, the fractional change in CE is very small as 
well. The aspect of classification in terms of outflow rate is be- 
ing looked into. The analysis is in progress and will be reported 
elsewhere. 
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